Dynamical SU(3) Linear a Model and 
the mixing of rj-r\ and a—fo mesons 



D. Kekez°, D. Klabucar 5 , M. D. Scadron c 

a Rudjer Boskovic Institute, P.O.B. 180, 10001 Zagreb, Croatia 

b Physics Department, Faculty of Science, University of Zagreb, Bijenicka c. 

Zagreb 10000, Croatia 

c Physics Department, University of Arizona, Tucson Az 85721 USA 



Abstract 

The SU(3) linear a model (LaM) is dynamically generated in loop-order using 
the nonstrange-strange basis. Only self-consistent logarithmic divergent graphs are 
needed, with quadratic divergent graphs replaced by SU(3) mass-shell equal splitting 
laws. The latter lead to an 7]'-r] mixing angle of 41.84° which is consistent with 
phenomenology. Finally this above SU(3)L<rM in turn predicts strong decay rates 
which are all compatible with data. 
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1 Introduction 



The original tree-level SU(2) spontaneous symmetry breaking (SSB) linear a model 
(LaM) interaction Lagrangian density is [1] 

L l ?} M = 9$ (O- + tl,T ■ 7T) V + ^(a 2 + 7T 2 ) - ^ 2 + 7T 2 ) 2 , (1) 



with tree-order chiral limiting (CL) couplings satisfying for f w ~ 93 MeV 

= ^, </ = ^ = A/., (2) 

Jit ^J-k 



here for quark fields, constituent quark mass m and cubic and quartic-meson couplings 
g', A. Although the g, g', A in (2) are not further specified in tree order, in loop order they 
are dynamically generated as 

g = -^L « 3.6276 , g = 2gm « 2.3 GeV , A = ^j- ps 26.3 . (3) 



We take m = 315 MeV, i.e., roughly one third of the nucleon mass Mjv, since this value 
most consistently satisfies the relations central for the present paper, namely mass shell 
equal splitting laws (MSESLs) considered below, in Sec. 3. (A dynamical quark mass 

- -i 1/3 

4 ™ s (-^)iGev ~ 320 MeV used in QCD [2] would lead to some 3% higher 



m = 
MSES 



3 

s value 4m 2 , whereas m = 325 MeV, used in Refs. [3], would lead to some 6% 
higher value.) While Refs. [3] recover the original chiral relations of Eqs. (2), the loop- 
order values in Eqs. (3) depend on relating slowly diverging log-divergent graphs with 
the more rapidly diverging quadratic-divergent graphs. 

In this paper we instead dynamically generate the SU(3) LaM using only the log- 
divergent graphs while replacing the quadratic-divergent graphs with the dynamical SU(3) 
MSESLs. In Sec. 2 we review the log-divergent gap equations in the CL and show their 
self-consistency in loop order. Then in Sec. 3 we replace quadratic-divergent mass gap 
equations with SU(2) and SU(3) MSESLs. In Sec. 4 we show how the SU(3) MSESLs 
lead to an 77-77' nonstrange-strange mixing angle of ~ 42°, in fact closely agreeing with 
the phenomenological value [4,5]. Finally in Sec. 5 we employ this SU(3) LaM with 
nonstrange-strange 77, 77' couplings to predict strong interaction decay rates for a^s ^n, 
a — > 7771-, / — > 7T7T and 77' — > rjirir, all in close agreement with data [6]. We give our 
conclusions in Sec. 6. 
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2 Self-consistent log divergent gap equations 



We begin with the non-perturbative loop-order equation for the pion decay constant 
fifn = fn in the soft-pion chiral limit [3]: 



1 = -^/ ' 



d A p 



(p 2 — m 2 ) 2 



(4) 



using the Goldberger-Treiman relation (GTR) m = f n g as in (2) with d A p = (2ir) A d A p 7 
where the quark mass m cancels out of this gap Eq. (4). This log-divergent gap equation 
(LDGE) (4) with g ~ 315 MeV/90 MeV ~ 3.5 requires an ultraviolet cutoff A « 750 MeV, 
separating the qq elementary particles n and m a ~ 650 MeV [6] with m nta < A from the 
bound-state qq mesons A < p(770), cj(780), ai(1260). This natural separation of LaM 
elementary particles from bound states is a consequence of the Z = compositeness 



condition [7] g 
Refs. [3]). 



In I \fW c or g = 3.6276 for N c = 3 (also dynamically generated in 



The self-consistency of loops "shrinking" to trees in the CL and their link to the LDGE 
are seen for quark triangle and quark box graphs. In the former case the quark triangle 
representing g' an7r is log-divergent with 



g' a ^ = -8ig 3 N c m J 



d A p 

(p 2 — m 2 ) 2 



2gm 



-AiN c g 2 



d A p 

(p 2 — m 2 ) 2 



2gm , (5) 



by virtue of the LDGE (4). Then using the quark-level GTR, Eq. (5) shrinks to the tree 
level g' anjT — > g' — m 2 T /2f 7T of (2) provided that m a = 2m, the Nambu-Jona-Lasino [8] 
(NJL) result also dynamically generated in Refs. [3]. Likewise the nn box graph in the 
CL gives the quartic quark coupling [3] 



A box = -8,iW/ *JL =2g 



(p 2 — m 2 ) 2 



-AiN c g< 



d A p 



(p 2 — m 2 ) 2 



= 2g\ 



(6) 



again via the LDGE (4). Then using the GTR, Eq. (6) becomes 



A 



box 



2g 2 



2gm = (/_ 

fir fir 



A 



tree 



(7) 



by virtue of the tree-level LaM couplings in Eq. (2). So again loops shrink to trees, while 
recovering the NJL scalar mass m a = 2m in this self-consistent fashion [3]. 

Next the CL quark bubble plus quark a tadpole graphs, although both being quadrat- 
ically divergent give a vanishing m 2 = in the CL (as required) provided the couplings 
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satisfy g' a ^ = m 2 r /2f. K , again recovering Eq. (2), but independent of the quadratically 
divergent scale. 

Lastly the pion quark triangle photon graph automatically normalizes the form factor 
F n (q 2 = 0) = 1 as expected. Specifically this quark triangle predicts [9] 

F n (q 2 ) = -iAN c g 2 [ dx [ — — , 

J J [p z — m l + x(l — x)q l \ l 

recovering F n {q 2 = 0) = 1 due to the LDGE (4). 



3 Mass-Shell Equal Splitting Laws 

To continue circumventing the dangerous quadratic divergent tadpole graphs, we first 
invoke the Lee null tadpole sum [10] characterizing the true (not false SSB) vacuum. 
Using only dimensional analysis, the vanishing tadpole sum requires [3] N c (2m) 4 = 3m*, 
or N c = 3 for the SU(2) LaM since we already know from Sec. 2 that the NJL relation 
m a = 2m is also valid in the LaM in the CL, as is g = 2ix / \f?> [7]. 

Away from the CL this NJL condition becomes for m « M N /3 ~ 315 MeV, 

ml - m 2 n = Am 2 ss 0.397 GeV 2 . (8) 

In Sec. 4 we will show that r/'-r] mixing requires a nonstrange-strange (NS-S) mixing angle 
(f) P w 41.84°, which in turn fixes the eta NS mass to be m NS = 757.9 MeV. Then the SU(3) 
extension of the mass-shell equal splitting law (MSESL) Eq. (8) is for m ao 984.8 MeV, 

< - m ms « °- 395 GeV2 . (9) 

which we again identify with the NS quark mass factor 4m 2 in Eq. (8). 

MSESL (9) follows from the empirical m 2 and from rn 2 s extracted in a phenomeno- 
logical way in Sec. 4 (see also Ref. [11]), but we have yet another way to avoid quadratic 
divergent amplitudes and evaluate the difference of the a and r] NS self-energies explicitly. 
It is encouraging that this explicit calculation below yields results which are reasonably 
close to MSESL (9), even though we consider only the lowest-order LaM self-energy 
graphs. Also, in counterdistinction to Eq. (9)Q, these graphs do not capture (at least not 

1 In Eq. (9), we plug in m^s = 757.9 MeV which does contain [11] the shift due to the gluon 
anomaly. 
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'H NS , , 'H AIS 



Fig. 1. Bubble graphs and tadpole graphs 

fully) the effect of the gluon anomaly which influences strongly the masses in the rj — rf 
complex. 

This other way to avoid quadratic divergent amplitudes is to subtract mass shell bubble 
(and tadpole) graphs of Fig. 1. Since the quadratic divergent tadpole graphs of Figs. 1 
b,d clearly cancel (due to the LaM coupling relation g' aoaQ(T = 9' msms a)i the remaining 
quadratically divergent bubble graph difference (Figs. 1 a,c) give the formal result 



m 



ao 



m VNS = 



_L 

-8iN c g 2 J dx J d 4 p 



p 2 — m 2 ao x{l — x) + m 2 
{p 2 + m^xil — x) — m 2 y 



x 



m 



(p 2 + m 2 x(l — x) — m 2 ) 2 



(10) 



Here we have combined propagator denominators using Feynman's trick 

1 j dx 
ab J [ax + 6(1 — x)] 2 



In other words, we evaluate the bubble graphs with propagator momentum p — > p — qx 
for mass shell values q 2 = m 2 ao , m 2 s for Figs. 1 a, c. Specifically for m ao = 984.8 MeV 
and m VNS = 757.9 MeV and a constituent nonstrange quark mass of 315 MeV, a computer 
calculation detailed in the Appendix, evaluates Eq. (10) as 

m 2 -m 2 = 5.83 m 2 -% 3.83 m 2 . (11) 
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We neglect the negative imaginary part of Eq. (11) (compatible with unitarity) since 
quarks in the bubble graph of Figs. 1 should be confined. Considering the delicate can- 
cellation due to the mean value (6 Jq x(l — x) — 1) giving in the second denominator 
in Eq. (10) m 2 NS x(l — x) — m 2 ss —0.035 m 2 , and especially recalling i) that we have 
used only the lowest order graphs and ii) that it is known that the gluon anomaly shifts 
m VNS upwards [which is taken into account in the result (11) only partially and indirectly, 
through the mass shell value q 2 = m 2 NS = (757.9 MeV) 2 ] so that the difference (11) is 
understandably somewhat overestimated, we suggest that the real part of (11) is not far 
from the numerical values of Eqs. (8) and (9), i.e. [12], 

m 2 aNS - m 2 = m 2 ao - m 2 ms = 4m 2 » 0.397 GeV 2 . (12) 

Stated another way, combining the partial fraction integrands in Eq. (10), we note that 
the leading p 6 terms in the numerator of (10) exactly cancel since they are quadratically 
divergent terms. Moreover, the log-divergent p 4 terms in (10) are proportional to 6(m 2 — 
m 2 NS )x(l — x) — m 2 ; they cancel using the mean value, resulting in m 2 — m 2 NS = 4m 2 . 
Thus again we support the MSESLs of Eq. (12). 



4 Dynamical i]'—r] Mixing 

Given the MSESL Eq. (12), one may extract the NS eta mass as m VNS ~ 757 MeV. 
Alternatively we may express the eta NS-S mass matrix as 



m VNS 


7 




m 2 





7 


m vs_ 


<t>P 





m 2 ,_ 



(13) 



where the NS-S pseudoscalar mixing angle (f)p determines the mixing relations 

\v) = cos0 P |?7 N s) - sin0 P |?7 S ) , \rf) = sm(f) P \r] NS ) + cos0p|r? s ) • (14) 

The angle 4>p is uniquely determined via the trace constraint 

m Ls + m vs = m l + m v> ~ L217 Gey2 > ( 15 ) 

(because the diagonal masses m v and are measured [6]) with off-diagonal hamiltonian 
matrix elements vanishing (r)'\H\rj) = (r)\H\r) f ) = 0, giving 

~m 2 NS = cos 2 0p m 2 + sin 2 0p m 2 , , m 2 s = sin 2 (f) P m 2 + cos 2 0p m 2 , . (16) 
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The two-level quantum mechanical solution of Eqs. (13), (15), (16) is the angle 



4>p = arctan 



m n-s m i =41 . 84 i 



rat 



ra. 



(17) 



Vns 



with masses 



m 



VNS 



757.9 MeV , m„ s = 801.5 MeV 



(18) 



Not only is the eta NS mass in (18) close to 757 MeV found via the MSESL Eq. (12), 
but the mixing angle in (17) is precisely the dynamical angle obtained via nonperturbative 
QCD. More specifically, Refs. [13,11] predict <pp as 



arctan 



(m 2 v , - 2m\ + m\){m\ - ml) = ^ 
(2ra 2 K — ml — m 2 )(m 2 , — ml) 



(19) 



found from the nonperturbative QCD gluon quark annihilation strength 

(m 2 , — m 2 )(m 2 — ml) „ 
(3 = y i \ — HL ~ 0.278 GeV 2 

4{m 2 K - ml) 



(20) 



and a constituent quark mass ratio X pa 0.78 ~ m/m s obtained from the NS-S QCD 
mass matrix [13,11] 



ml + 2(3 



V2[3X 




ra 2 


-ml + (3X 2 


<pp 


m 2 , 



(21) 



Equations (18), (20), (21) have the solution 



tan(20 P ) = 2V2(3X{ 



in 



?IS 



m 



VNS 



y 1 = 9.02 or P = 41.84° . 



(22) 



Again we see that 4>p in (17), (19) and (22) are extremely close in magnitude. Moreover, 
the dynamical approach to 77-77' using the Schwinger-Dyson (SD) and Bethe-Salpeter (BS) 
integral equations found (pp ~ 42° (that is, in terms of the singlet-octet state mixing 
angle 6 P = <f) P — arctan \[2 pa —12.7°) [14]. Its subsequent refinement [11] also included 
the effect of the "strangeness attenuation parameter" X in the SD-BS mass matrix. The 
SD-BS estimate was X = 0.663, again close to the constituent quark mass ratio m/m s 
found there. Fitting the trace constraint (15) then led to (3 — 0.277 GeV 2 , practically the 
same as Eq. (20), to ra VNS = 757.87 MeV and m vs = 801.45 MeV, almost the same as Eq. 
(18), and to P = -13.4°, that is, <p P = 41.3°, very close to <f> P in Eq. (22). (These results 
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were for the original parameters of Ref. [14]. Reference [11] also varied the parameters to 
check the sensitivity on SD-BS modeling, but the results changed little.) Because of the 
close link between Eqs. (17), (18), the QCD Eqs. (19), (22) and the SD-BS scheme, we 
suggest that (ftp ~ 41.84° is the dynamical r]'-r] mixing angle in the NS-S quark basis. It 
corresponds to [4,5] the singlet-octet angle Op = (ftp — arctan \[2 ps —12.9°. Also note that 
Ref. [11] showed there is no contradiction between our approach utilizing one state-mixing 
angle, and the mixing scheme employing two angles pertaining to the mixing of the decay 
constants (see Refs. [5], esp. the second reference for review.) Not only is the difference 
small in the NS-S basis, but our Ref [11] also showed that our results are in agreement 
with what is found in the two-angle scheme [5]. 

It is also satisfying that the phenomenological analysis of the NS-S rj'-r) NS mixing 
angle extracts [4] (ft P = 43.2° ±2.8° from T — > PP decays, (ft P = 36.6° ±1.4° from V -> P7 
and P -> V'y decays, <p P = 41.3° ± 1.3° from P -> 77 decays and <p P = 40.2° ± 2.8° from 
J/ift — > pi], pi]' and urn decays. Moreover the recent Refs. [5] obtain (ftp = 39.3° ± 1.0° by 
global phenomenological fits and (ftp = 42.4° as their theoretical prediction, which is all 
within the region of the dynamical (ftp angles in Eqs. (17) or (19). 



5 SU(3) LctM Strong Decay Rates 



We have thus far used the LDGE (4), induced the MSESLs m 2 NS — m\ = m 2 ao — m 2 NS = 
m 2 —m 2 K = Am 2 ss 0.397 GeV 2 (for k(805— 820) advocated by, e.g., Delbourgo and Scadron 
[12]) and the NS-S mixing angle (ftp ~ 41.84° all while avoiding quadratic divergent 
graphs and extending the SU(2) LaM to SU(3). In the latter case the cubic meson LaM 
Lagrangian density has the SU(3) form [12] 

L^ c = d ijk (g' S PpS*P>P k + gfsss&SlS") . (23) 



Then with f n ^ 93 MeV and m » Mjv/3 ^ 315 MeV, the MSESLs above suggest the 
Lagrangian g' SPP couplings 



2 2 

Tfl — Tfi 

4m = - 2.13 GeV , (24) 



2 2 
mi„ — mi 



2U 

m 2 

9'a„= a ° y^ 3 « 2.13 GeV , (25) 
9' k Ktv — — ^y— ^ ~ 2.13 GeV , (26) 



along with g' aoV7T = cos(ft P g' aoms7T , g f rfao7T = sm(ft P g' aoms7T , etc. 
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The nonstrange a decay rate is predicted as [6,15] 

T(a NS -> tttt) = ^(2^ s7r7r ) 2 g-|L w 754 MeV , (27) 



for m aNS 650 MeV and \p\ = 294 MeV. This rate is compatible with Weinberg's mended 
chiral symmetry estimate [16]: 

T^s « « 676 MeV . (28) 



Likewise the SU(3) LaM ao — > ^ decay rate is 

\P\ r_ . 



^LaM{.CLQ -> 7/71") 



87rm2 Q 



[2^ 



a-OVNS 



n COS </> P 



133 MeV 



(29) 



for p = 321 MeV, g' aoms - K ~ 2.13 GeV, 0p = 41.84°. One may infer a nearby ao rate from 
the PDG tables [6]. Specifically the rate ratio 

T ^^ KK } = 0.177 ±0.024 (30) 
r(a ^7r) V ; 



and T(a -> M) « 24.5 MeV from Refs. [17], then suggests T(a -> ^vr) « 138 MeV, 
near Eq. (29). Also, this predicted LaM decay rate (29) is not too distant from the high 
statistics data [18] 

T aovn = (95 ± 14) MeV . (31) 

The SU(3) companion / (980) — > rnr rate is estimated [6] to be 

T(f inr) « (47 MeV) (0.781) « 37 MeV (32) 



assuming the small r(/ 77) ~ 0.56 keV rate in the 1998, 1996 PDG tables combined with 
the measured branching ratio £>( i /o77) ~ 119 x 10~ 5 . On the other hand we must account 
for scalar a-fo mixing (the analogue of pseudoscalar r}-rf mixing). Thus in the NS-S basis 
we define in parallel with Eq. (14) 

\a) = cos0 5 |o- NS ) - sin0 5 |o- s ) , |/ ) = sin0 5 |o- NS ) + cos(f) S \a s ) , (33) 



and estimate 0s from the measured decay rate ratio 

r(/ovr7r) 3 /470 MeV\ /sin0 5 \ 2 37 MeV ,, . n , x 

-7^ — '- « - — « « 0.39 or 6 S « 18.3° . 34 

r(a ?77r) 2 V321 MeVMcos0 P y 95 MeV V ; 
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Prior theoretical estimates were \(j> a \ ~ 16°, 20° [12] and 14° [19]. The DM2 data of 1989 
[6] also suggests from J/ip — ■> utttt that /o(980) is mostly ss (not nonstrange), compatible 
with (34) (and near the <fi (1020) which is known to be almost all ss) [12,20]. 

Lastly we calculate the strong decay rate rj — > rj7nr in the context of the SU(3) 
LaM [21], with a ,a, f poles contributing as rj' — > a 7r — > rjim (4 modes), rj' — > r\a — > 
r]7i7i,i]' — > 77/0 — > r77T7r. Although the 4a pole modes should dominate, the well-known 
LcxM rj — > ?77r7r contact term 3A [normalized to the quartic term in the SU(2) Lagrangian 
Eq. (1)] has the opposite sign relative to do, o and fo poles and "miraculously cancels" 
them [22] due to chiral symmetry - assuming one treats the a , a, / poles in narrow 
width approximation. While T ao /m ao , Tf /mf ~ 1/10 as needed, the a is broad with 
Y a /m c ~ 1. 

Then after the chiral cancellation, we must still account for the broad-width a inverse 
propagator as s — m 2 a + im a T a with \s — m 2 a \ « ^777,0-^1. Thus the net rj' — > rjn 7r LaM 
amplitude has the magnitude 

WTUv' - WV)| « |^EL| „ |^| „ 5.7 . (35) 

Here we [11] estimated g'^ m cos^p sin^p^^. m 1.06 GeV. Then the net SU(3) LaM 
decay rate is predicted to be (folding in the 3-body phase space integral [23]) 

IW?/ -> ^7T°7r ) = 1.06|M-V| 2 keV » 34.4 keV . (36) 

A slight increase of this rate (36) is due to the 10% non-narrow widths of the ao and fo 
poles. Recent data gives [6] T(r}' — > ^7r 7r ) = (42 ±4) keV. The total decay rate assuming 
isospin invariance is 

^LaAlW ? ? 7r7r ) = F L aM(v' ~^ ^71°) + T LaM (rj' -> rjH~H + ) 

= 3 x (34.4 ± 4) keV = (103 ± 12) keV , (37) 

near the total observed rate of 3x (42±4) = (126±12 )keV. We know of no other dynamical 
scheme (such as using the original singlet-octet mixing angles [21]) which recover all the 
approximately needed SU(3) strong decay rates (27), (29), (36), (37) as found above. 



6 Conclusion 

In this paper we have consistently avoided dealing with quadratic divergent graphs when 
computing SU(2) and SU(3) linear a model {LaM) diagrams. Instead in Sees. 2 and 
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3 we work only with self-consistent log-divergent gap equation integrals Eqs. (4), (10). 
Sections 3 and 4 extend this pattern from SU(2) to SU(3) dynamical mass-shell equal 
splitting laws, leading to the off-diagonal eta nonstrange and strange constituent quark 
masses m VNS ~ 757.9 MeV and m vs ~ 801.5 MeV. Then the dynamical rj'-r] mixing angle 
in the NS-S basis is (f>p ~ 41.84° compatible with nonperturbative QCD and near many 
phenomenological analysis of this NS-S angle (see, e.g., Refs. [4,5,13]). 

Stated another way, the only SU(3)-breaking pattern we allow is characterized by the 
constituent quark mass GTR ratio [24] as used in the phenomenological analysis of Refs. [4] 
m s /m ss 2f K /f n - 1 S3 1.44 for f K / U ~ 1-22 as measured [6]. Then in Sec. 5 the SU(3) 
SPP LoM couplings (again following the above MSESLs) of Eqs. (8), (9), (11), (12) in 
turn predict strong interaction a^s 7l71 , a o fo 71-71 v' ~^ V 7171 decay rates all 

compatible with data [6]. 
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Appendix: On the bubble graph integral 

If the integrand of Eq. (10) is rewritten using the common denominator D(x,p 2 ), as 
f(x,p 2 ) = N(x,p 2 )/D(x,p 2 ), one should note that the 0(p 6 ) terms in its numerator 
N(x,p 2 ) cancel exactly. The numerator is thus a polynomial of degree 2 in (p 2 ): N(x,p 2 ) = 
Cq(x) + Ci(x)p 2 + C2(x) {p 2 ) 2 ■ The integrand is therefore conveniently written as the sum 

f(x,p 2 ) = £ f, t (x, p 2 ) = j: Cl (x) JfL^ . (38) 

The four-dimensional integral over p is effectively one-dimensional because the integrand 
depends on p 2 only. After the Wick rotation, we performed this integration analytically, us- 
ing the Mathematica program package. The log-divergent integral / d 4 p f2(x,p 2 ) depends 
on our ultraviolet cutoff A = 750 MeV required by Eq. (4). After the p 2 -integration, the 
logarithmic forms 

, I m 2 — m 2 (1 — x)x\ , 

l(x) = In — ^ ±- 39 

\m z - m^ o (l — x)x J 

appear in the integrand, requiring some care. The mild divergences at the points Xo = 
0.115698, xi = 0.222046, x 2 = 0.777954, and x 3 = 0.884302 correspond to the roots of 
polynomials x i— > m 2 — m 2 s (l — x)x and x i— > m 2 — m 2 (l — x)x. In order to perform 
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the residual x integration of the functions x \— > / d A p fi(x,p 2 ) (i = 0, 1,2), the interval 
[0,1] is divided into five integration regions, [0, x±], [xi,X2], [^2,^3], [^3,^4], and [X4, 1]. 
These integrations were numerical, with an adaptive algorithm which can handle the mild, 
integrable singularities appearing at the edges of the integration regions. 
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